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Abstract 

The generalized Bretherton equation is studied. The classifica- 
tion of the meromorphic traveling wave solutions for this equation is 
presented. All possible exact solutions of the generalized Brethenton 
equation are given. 

1 Introduction 

At present there exists a lot of methods for finding exact elliptic solutions 
of autonomous nonlinear ordinary differential equations. Let us name only 
a few: the Weierstrass function method [I1I2], the Jacobi elliptic-function 
method [3HS] and their different extensions and modifications [BHI2] ■ Making 
use of such a method one can not but come across the following questions. 
Whether all families of elliptic solutions are found. Whether an equation 
does not have elliptic solutions at all, if a method fails to find any. This 
questions are addressed very seldom (nevertheless see [T3HI5])- 

The aim of this article is to present an algorithm, which enables one to 
find all families of doubly periodic meromorphic solutions satisfying an au- 
tonomous nonlinear ordinary differential equation. We will use an approach 
suggested in [TH|T5] . As an example we take the following third order differ- 
ential equation 

Wzzz + ^V^Wzz + 11^1 + Pwzz + 2aw'^ + 5wWz + "^Wz + en? + iiw + = 

(1) 

with 77 7^ 0, a 7^ 0. We find all conditions for elliptic solutions to exist and 
construct elliptic solutions in exphcit form. 
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This article is organized as follows. In section 2 we present our method 
and give explicit expressions for elliptic solutions of an autonomous nonlinear 
ordinary differential equation. In section 3 we consider an example and 
classify elliptic solutions of equation ([1]). 



2 Method applied 

Consider an algebraic autonomous nonlinear ordinary differential equation 

E[w{z)] = 0. (2) 

Let us look for its elliptic solutions. If w{z) is such a solution, then equation 
([2]) has the family of elliptic solutions w{z — Zq) with Zq being an arbitrary 
constant. Equation ([2]) necessary possesses an elhptic solution, if it admits at 
least one Laurent expansion in a neighborhood of the pole z = zq. Without 
loss of generality we may build Laurent series in a neighborhood of the point 
z = Q 

p oo 

^(^) = j2^ + Y.'"^^'^ 0<|z|<£. (3) 

k=l k=0 

Here p is the order of the pole z = 0. 

Proposition. Suppose Laurent series ([3]) with uniquely determined co- 
efficients satisfies equation then this equation admits at most one mero- 
morphic solution having a pole z = with Laurent series (E]) . 

This proposition follows from the propertied of Laurent series and unique- 
ness of analytic continuation. As a consequence, equation (|2]) may have at 
most one elliptic solution possessing the pole z = with Laurent series ([3]). 
The order of an elliptic function is defined as the number of poles in a par- 
allelogram of periods, counting multiplicity. 

Our algorithm for finding elliptic solutions of equation ([2]) in closed form 
is the following. Note that we omit arbitrary constant zq. 

Step 1. Perform local singularity analysis around movable singular points 
for solutions of equation ([2]). 

Step 2. Select the order M of w{z) and take K distinct Laurent series 

P^ (i) OO 

^^^n^) = E^ + E4^^^'> 0<|z|<£„ t = l,...,K. (4) 

k=l k=0 

from those, found at step 1, in such a way that the following conditions 

Ec«=0, f2p^ = M. (5) 
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hold. 

Step 3. Define the general expression for the elliptic solution w{z) pos- 
sessing K poles ai, . . ., Oii- in a parallelogram of periods (see [HHTG]). The 
Laurent expansion in a neighborhood of the point z = ai is w^^\z — ai). In 
other words, take the following expression for w{z) 




Here p{z) is the Weierstrass function satisfying the equation 

{p, f = 4p^ - g2P - 93, (7) 

C{z) is the Weierstrass C^function, Hq is a constant. 

Step 4- Find the Laurent series for w{z) given by (j6]) around its poles 
Oi, . . ., aK- Without loss of generality set ai = 0. Introduces notation 

Ai =^ p{ai), Bi =^ p^(ai), i = 2, K, and 

K 



ho = ho-^ c^l\C{ai) - ^ c^llp{ai). 



1=2 iG / 



With the help of addition formulae for the functions p and ( (see 
rewrite expression ([6]) as 



wiz] 



(A;-l)! 



^ cg(p.(.)+i?.) 

2(p(^)-A,) + (fc-1)! rfz^ 




(9) 



iStep 5. Compare coefficients of the series found at the second and the 
fourth steps. Form a system of algebraic equations. Add to this system 
equations Bf = AAf — g2Ai — g^, i = 2, . . ., K. The number of equations 
in the system is slightly more than the number of parameters of elliptic 
solution ([9]) and equation (|2]). Solve the algebraic system for the parameters 
of the elliptic solution w{z), i.e. find ho, g2, gs, Ai, Bi, i = 2, . . ., K. In 
addition correlations for the parameters of equation ([2]) may arise. If this 
system is inconsistent, then equation does not possess elliptic solutions 
with supposed Laurent expansions around poles. 

Step 6. Check-up obtained solutions, substituting them into original 
equation. 
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In the case (^f — 27 gl = the elhptic function p{z) degenerates and con- 
sequently elhptic solution degenerates. The invariants g2, gs are related 
with the half-period ui, U2 by means of equalities 

_ 60 _ >^ 140 

(n,m)7^(0,0)^ ^ (n, m)^^ (0, 0) ^ ^ ''^ 

With the help of presented algorithm one can construct any elliptic so- 
lution of equation ([2]). Note that if equation ([2]) possesses only distinct 
Laurent series in a neighborhood of poles, then the orders of its elliptic so- 
lutions is not more than ^^iPi, where pi {i = 1, . . ., N) are the orders of 
poles y^|15j. Thus we see that our approach may be used, if one need to 
classify families of elliptic solutions satisfying equation ([2]). 



3 Elliptic solutions of an autonomous third 
order ordinary differential equation 

As an example let us classify doubly periodic meromorphic solutions of the 
third order nonlinear ordinary differential equation ([1]). Without loss of gen- 
erality we may set r/ = 1, cr = to obtain 

Wzzz + SwtOzz + wl + jSwzz + 2aw^ + 6wWz + + + z/ = 0, (11) 

where a 7^ 0. Our results are gathered in theorem O 

Theorem 1. Equation fflTl) admits exact elliptic solutions if and only if 
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Case I : 5= —a, 7 = -^"/? (12) 



54^3/3^ + 783a^l3^n + 3222a/3V^ + 2816^^ + lOSOOOai^^ 
-4725a2^V - 32ma(3ufi ^ 



30^ _ 2a(990a + 833/?) f (^"^ , 15480720 2 



Case II : S = a, 7 = , n = —2a \ a 

77 ' ' 3773 ' ^ V 3 3195731 

2916 \ /4/33 88719651 26053515945 o 4005 ^2 

+— — a/3 , 1/ = -a + a^(3 + a"^ + —— a/S'^ 



539 y ' V 27 12782924 ' 6889996036 1232 

(14) 

5929 .5929 

^ 5634 ^ 103230 ^' 
104796 ^2 5558014908 ^ 2 21448672012224 . , ^ ' 

^ 4235 ^ 125546575 ^ 1042111900445 
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Proof. Equation flTT]) possesses two different asymptotic expansions cor- 
responding to Laurent series in a neighborhood of poles. They are the fol- 
lowing 



, 11 115 -12a 324^2 _ 309a(5 + 11(52 + 1083a/? 

w^^>{z) = + — — + — — + . . . 

^ ^ az^ 19az 8664a 

(16) 



yji^)(z) = -^+ (l6-^a-l(3] + ... (17) 



7z \7 49 3 
All not written out coefficients of these series are uniquely determined. With- 
out loss of generality we construct Laurent series in a neighborhood of the 
point z = 0. Thus we see that equation ffTTjl may have an elliptic solution of 
order two and an elliptic solution of order three. The second order elliptic 
solution possesses one pole inside a parallelogram of periods. The Laurent 
expansion in a neighborhood of the pole z = is w^^\z). While the third 
order elliptic solution has two distinct poles inside a parallelogram of periods. 
The Laurent expansion in a neighborhood of the pole 2 = is w^^\z) and 
the Laurent expansion in a neighborhood of the pole ^ = a is w^'^\z — a). 
All other elliptic solutions are obtained with a help of transformation z 1— )■ 
z — zq, where Zq is an arbitrary constant. Necessary condition for the second 
order elliptic solution to exist is cL^] = and we obtain the first constrained 



Wiz) = ^r- + ho -\ ^^Z H ^7^2 + \^'^) 



between the parameters in expression f[T2|) . Necessary condition for the third 
order elliptic solution to exist is clI + clI = 0. Consequently we get the first 
constrained between the parameters in expression (fT4|) . The second order 
elliptic solution can be written as 

w{z) = c^11p{z; g2,g3) + ho. (18) 
Expanding this function in a neighborhood of the pole z = 0, we obtain 

^ ' ■ 20 ' ' 28 
Comparing coefficients of the series (|T6l) and (|T9|) . we find the parameters ho, 
g2, Qz of elliptic solution ( IT8|) 

15 a(32/x + 3/32a) a^ (320z/ - 7/?=^a - 48/?u) 

8 880 38720 

and the second constraint in f|T2|) . In fact we take seven coefficients of each 
series. Excluding the case of degeneracy, i.e. §2 — 27gl = 0, we force the 
condition f|T3|) . The third order elliptic solution is the following 

wiz) = c%iz) + c«(C(^) - C(^ - a)) + ho- c«C(a), c^'l = -c^-l (21) 
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where p{z) = p{z; g2, g^), C{z) = ({z]g2,g3). Expanding this function in a 

def 

neighborhood of the poles z = 0, z = a and introducing notation A = p(a), 
B =^ pz{a), we get 

t \ ^-2 C_i (1 I C_2fi'2 C_i^ I 2 , 

^(^) " ^ + ~ + + + ] ^0 2~ ^ ^ ^ ^ 



(2) 

«^(.) = ^ + /^o + c^A + {c^lB - c^\a} e + {c^l {SA^ - f ) 



(23) 



Comparing coefficients of the series with coefficients of w''^\z) and 

— a) accordingly, we obtain an algebraic system. In this case we take 
five coefficients of each series. Solving this system together with equation 
B^ = AA^ — g2A — g3, we find the parameters ho, g2, g3, A, B 

111a P 2 f 8139 ^ 43244955 . ^^x 

ho = 1 — , On = a pa H a -\ 

o.r. y2 \^ 23716 ^ 140612164 48 J 



2156 



3 , 3163 2 15684093 2 19911334005 3 /3'^ 



189728' 562448656' 1667379040712 1728 



/ /3 69a \ „ 2 

A = a\ — ^ , 5 = -or 

V24 2156/ V 



15/? 43119a\ 
308 ^ 913066 J 



(24) 



and the second and the third correlations in (1141) . Rewriting elliptic solution 
fl2T]) in terms of the parameters A, B, yields 

w{z) = c_2p{z) - _ ^) + ho^ (25) 

In order to exclude the case of degeneracy we should force the condition flTBl) . 

Substituting obtained elliptic solutions into equation fITT]) . we see that 
this equation indeed possesses solutions of the form ( IT8|) . (I25|) provided that 
the parameters a, /3, 7, 5, /i, v satisfy correlations (fT2|) . ( 1T4|) accordingly. 

Note that if conditions ( !T3|) . ( [T5|) are not valid, then expressions (fT8l) . 
(1^ give simply-periodic or rational solutions of equation (ITT]) provided that 
correlations (11 2p . (1141) hold accordingly. 

In conclusion we would like to note that solutions of equation ([1]) with 
cT = provide traveling wave solutions w(x, t) = h + w{z), z = x — {-y — 6b)t 
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of the following third order nonlinear partial differential equation 

Ut = Uxxx + ^V^Uxx + '7^x + '^2Uxx + 2aM^ + 5uUx + aiu^ + a^u, (26) 

where 02 = /? — 36, ai = — 6a6, = 6a6^ + /i and & is a root of the equation 
2ah^ + fib — u = 0. Along with this, solutions of equation ([T]) give traveling 
wave solutions u{x,t) = w{z), z = x — fj,t of the fourth order nonlinear partial 
differential equation 

ut = dx {uxxx + ^VUUxx + V^l + f^'^xx + 2au^ + SuUx + ^Ux + au^] . (27) 

In this case the parameter v in equation ([1]) appears as a constant of inte- 
gration. 

4 Conclusion 

In this article we have studied the problem of constructing exact elliptic solu- 
tions of autonomous nonlinear ordinary differential equations. Our approach 
is based on comparing coefficients of Laurent expansions found by asymp- 
totic methods and coefficients of Laurent series for the general expression of 
an elliptic function. The method presented in this article generalizes sev- 
eral other methods with an a priori fixed expression for an elliptic solution. 
Our approach may be applied if one needs to classify elliptic solutions of an 
autonomous nonlinear ordinary differential equation. 
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